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Abstract 

We explore the possible impact of galactic and ex- 
tragalactic foregrounds on measurements of the cosmic 
microwave background (CMB). We find that, given our 
present understanding of the foregrounds, they are un- 
likely to qualitatively affect the ability of the MAP and 
Planck satellites to determine the angular power spec- 
trum of the CMB, the key statistic for constraining cos- 
mological parameters. Sufficiently far from the galac- 
tic plane, the only foregrounds that will affect power 
spectrum determination with any significance are the ex- 
tragalactic ones. For MAP we find the most trouble- 
some foregrounds are radio point sources and the ther- 
mal Sunyaev-Zeldovich (SZ) effect. For Planck they are 
these same radio point sources and the Far Infrared Back- 
ground. Prior knowledge of the statistics of the SZ com- 
ponent (either via theoretical calculation, or higher fre- 
quency observations of just a few percent of the sky, such 
as will be done by balloon-borne experiments) may signif- 
icantly improve MAP's determination of the CMB power 
spectrum. We also explore the foreground impact on 
MAP and Planck polarization power spectrum measure- 
ments. 

Subject headings: Cosmology, Cosmic Microwave Background 
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Introduction 



Much attention has been paid in recent years to the 
nature of the foregrounds which obscure our view of the 
background. This attention has resulted in discoveries 
about the nature of the foregrounds as well as meth- 
ods for estimating CMB anisotropy from foreground- 
contaminated data. From studying these developments, 
we have concluded that for planned large area, multi- 
frequency experiments, such as MAPQ and Planck^, the 
foregrounds are unlikely to be responsible for qualitative 
degradation of the primary cosmological results. 
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This happy situation is due to a number of factors. 
First, there is a window in frequency space, where, at 
high galactic latitude, CMB fluctuations are the brightest 
diffuse source in the sky. Second, the high-latitude galac- 
tic foregrounds are very smooth; they do not have much 
small-scale fluctuation power. Third, foregrounds, un- 
like instrument noise and some systematic error sources, 
are suppressed at small angular scales by the finite an- 
gular resolution of the telescope. Fourth, point source 
count estimates suggest that only a small fraction of pix- 
els in the MAP and Planck maps will be affected — and 
these can be identified with threshold cuts and removed. 
Finally, even if uncertainty in a particular mode of the 
CMB map is dramatically increased by the presence of 
foregrounds, uncertainty in the CMB power spectrum 
may not be significantly affected. This is due to the 
fact that, at least in the foreground free case, the dom- 
inant source of power spectrum uncertainty (except at 
the smallest angular scales) comes from sample variance, 
not instrument noise. 

The primary cosmological results — determination of 
cosmological parameters — depend mostly on how well 
the power spectrum is measured. We thus focus on 
the impact of foregrounds on the determination of this 
power spectrum. Our method for estimating this im- 
pact can be considered to be a generalization of those 
based on Wiener filtering by Bouchet, Gispert, and Puget 
(1995, hereafter "BGP95") and Tegmark and Efstathiou 
(1996, hereafter "TE96") as weU as that of White (1998, 
hereafter "W98"). All these approaches take the CMB 
and foregrounds to be statistically isotropic, Gaussian- 
distributed fields. Given this assumption, estimation 
of the power spectrum errors is straightforward, as de- 
scribed below. 

As is always the case with parameter estimation, how 
well the desired parameters can be recontsructed depends 
on the assumed prior information. The methods of TE96 
and W98 essentially assume that the foreground power 
spectra are known with infinite precision a priori. The 
most important difference between our method and theirs 
is that we only assume finite prior information about the 
foreground power spectra. 

Although the method of BGP95 was derived assuming 
Gaussianity, they have tested it with non-Guassian sim- 
ulations of Planck Surveyor maps (see, e.g., Gispert and 
Bouchet, 1996). Their results lend credibility to the fore- 
casts derived analytically under the Gaussian assump- 
tion. 

Below we first describe our methods for estimating 
the power spectrum uncertainties given an experimental 
configuration and foreground model. In section III we 
describe our model of the foregrounds, which is based on 
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that detailed in the Planck Phase A proposal (Bersanelli 
et al. 1996), and the High Frequency Instrument (HFI) 
and Low Frequency Instrument (LFI) proposals^. 

To date, foregrounds have been essentially ignored in 
estimates of the cosmological parameter uncertainties^. 
We find that they are unlikely to qualitatively change 
the results. Although for MAP this conclusion depends 
somewhat on the amplitude of the contribution from the 
Sunyaev-Zeldovich effect, which is not yet sufficiently 
well-determined. 

Not only does the amplitude of the SZ power spec- 
trum affect the ability of MAP data to constrain the 
CMB power spectrum, but so does our prior knowledge 
of it. This is fortunate, because while the amplitude is 
completely out of our control, we can do something about 
how well we know it. We emphasize that the prior infor- 
mation we need is not of the actual SZ map, but of the 
statistics of the map. The statistics can be calculated 
theoretically, or by actually measureing the SZ map over 
only a few per cent of the sky. 

Higher order moments of the probability distribution 
may also be of interest if the CMB statistics are non- 
Gaussian, which they will be to some degree even if the 
primordial fluctuations are Gaussian. Therefore, we also 
estimate how well the amplitudes of individual spherical 
harmonics can be determined. The uncertainty on these 
amplitudes is much more strongly affected by the pres- 
ence of foregrounds than is the uncertainty on the power 
spectrum. 

Effects due to contributions not in ones model of the 
data may be detrimental and our formalism does not 
address such a problem. Nevertheless, we find it encour- 
aging that for the quite general model we have chosen, 
where the data are required to simultaneously constrain 
thousands of foreground parameters, the results look very 
good. 

2. Methodology 

We assume that the experiment measures the full 
sky in = l,...,nch channels, and model the (beam- 
deconvolved, spherical-harmonic transformed) map data 
as due to the CMB, foregrounds and noise: 



(1) 



^These proposals are not yet publically available. Most of the work 
referred to here will soon be available as Bouchet and Gispert 
(1998). 

notable exception is the "LDB" case in Bond, Efstathiou & 
Tegmark (1997) which was based on calculations by the TopHat 
group of CMB pixel errors after pixel-by-pixel subtraction of fore- 
grounds in their Antarctic maps. 



where i runs over the components, (i = is CMB, i > 
are the foregrounds) and Qii, gives their frequency depen- 
dence. In the following we usually suppress all indices 
and use a notation in which Eq. |l| becomes: 

A = g^a + n. (2) 

Throughout we assume that the noise is spatially uni- 
form, Gaussian-distributed, and uncorrelated from chan- 
nel to channel. Therefore, W =< nn^ is given by 



(3) 



where the weight-per-solid angle for map v, Wi,, equals 
I / i'^t^pi^) ^ is the standard error in the temper- 
ature determination of a map pixel with solid angle 
ripix, and Bu i is the beam profile — which for a Gaus- 
sian beam with full-width at half-max V8 In 2cr IS given 
by exp{—{l<jY/2). The beam-damping of the weight ma- 
trix is due to the fact we are describing the noise in the 
beam- deconvolved maps. 

If we make specific assumptions about the statistics 
of the CMB and foregrounds then we can determine how 
well we can measure the parameters of those statistical 
distributions. For simplicity and specificity we assume 
the CMB and foregrounds to all be statistically isotropic 
and Gaussian-distributed. In this case a complete sta- 
tistical description of the data is given by the two-point 
function: 



M = (AAt) 

= g^(aat)g + W-^ 



(4) 



If, in addition to statistical isotropy, we assume that each 
of the foreground components are uncorrelated then we 
can write 

{^Im'A.'m') = CilSlliSmm'^ii' (5) 

and Eq. ^ simplifies to (with indices restored): 
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Given the data, we could write down and calculate 
the posterior probability distribution of the parameters, 
Cii, or any other parameterization, Op, of M. The poste- 
rior is proportional to the product of the likelihood and 
the prior. In the limit that the posterior distribution of 
Up is Gaussian, the expectation value for the covariance 
matrix of the parameters is given by the inverse of the 
"posterior" Fisher matrix. 
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Note that the trace is a sum over fe, ms and vs. M is 
block-diagonal with block size rich by rich so its inversion 
is readily feasible. The matrix F, or rather its inverse, 
is exactly what we want, the expectation value of the 
covariance matrix of the parameters. We are interested in 
calculating this parameter covariance matrix for various 
parameter choices - in particular the Cu - as well as 
assumptions about their prior distributions. 

We parameterize the (diagonal) prior as zero for i — 
and 

FSj' = HCu? (8) 

for i > where Cu are the assumed actual power spec- 
tra, to be discussed in the next section. Note that 
if we take the foreground CuS to be a priori perfectly 
known {a — s- oo), then Eq. |^ gives the Fisher matrix for 
the Wiener filter method of foreground removal (TE96, 
BGP95), an explicit expression for which is in W98. In 
the absence of foregrounds it is equivalent to that in 
Bond, Efstathiou & Tegmark (1997, hereafter "BET97") 
and for a single channel experiment it is equivalent to 
that given by Knox 1995 and by Jungman et al. (1996). 
Below we vary a to see quantitatively how the strength of 
our prior assumptions determines the ability to measure 
Coi. 

It is straightforward to generalize the above to include 
polarization information. Maps of the Q and U Stokes 
parameters can be decomposed into two components, 
and (Kamionkowski et al. 1997; Zaldarriaga & Sel- 
jak 1997), which are now in addition to the temperature 
component a^. In general, we can write the contribu- 
tion from each component as a^;^ and the data in each 
channel as A'^ where the superscript is either T, E 
or B. Then the covariance matrix for the data (Eq. ^) 
becomes 



error in the minimum variance estimate of a is 



M, 



lm,l'7n' 



—rSbb'S^u' 



Sll'Sri 



(9) 



where C^l' equals C^^ 
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for b = b' = T, 



b = b' = B, and = 
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um^zhn / for 6 = T, b' = E. 
All other elements vanish. Thus, while the matrix of 
Eq. ^ is block-diagonal in blocks of dimension rich, this 
matrix is block-diagonal in blocks of dimension Zrich- 
This approach generalizese the multi-frequency Wiener 
filter error forecasting of Bouchet et al. (1998, hereafter 
"BPS"), who generalized the single- frequency, no fore- 
ground, treatment of Zaldarriaga et al. (1997). 

We may also be interested in how well an individual 
mode can be measured. The covariance matrix for the 
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where we have assumed a prior probability for a that 
is Gaussian-distributed with weight matrix WP'''°''. For 
example, we may wish to assume that foreground i 
has variance GaSwSmm' in which case VF^^^,;,^, = 
^/CiiSiiiSfnm'- With this prior, this is the variance 
given by the Wiener filter procedure. Without the prior 
it is the variance given by the pixel-by-pixel subtrac- 
tion procedure of Dodelson (1996) and also of Brandt 
et al. (1994) (except for their non-linear parameter de- 
pendences). When there are more foregrounds than 
channels, g^Wg is singular and therefore addition of 
a prior is necessary to make {SaSa^) finite. For more 
flexibility in the prior choice later, we define /? so that 
^iiniiim — P/^ii ■ Notc that Eq. |l^ docs not assume any- 
thing about the statistical properties of the foregrounds 
and CMB — except through the prior, which we have ex- 
plicitly assumed to be Gaussian. 

3. Foreground Models 

Our foreground model is based on that developed for 
the Planck Phase A proposal (Bersanelli et al. 1996) 
and updated in the HFI and LFI instrument proposals. 
We refer the interested reader to these proposals and to 
Bouchet and Gispert (1998, hereafter "BG98"). Below 
we briefly describe our model, with an emphasis on the 
modifications and additions we have made. In all cases, 
these alterations make the model more pessimistic. 

3.1. Galactic 

Analyses of the DIRBE (Diffuse Infrared Background 
Explorer) and IRAS (Infrared Astronomy Satellite) Sky 
Survey Atlas maps have determined the shape of the dust 
power spectrum to be C; ex l~^ (Gauttier et al. 1992; 
Wright 1997) or Ci cx r^-^ (Schlegel et al. 1998). We 
assume C/ cx l~^'^ since it is the more pessimistic choice, 
given that we normalize at large angles. 

We take the same C/ shape for the free-free power 
spectrum because both the dust intensity and free-free 
are expected to be from the same warm interstellar 
medium. Indeed, there is strong observational evidence 
for a correlation (Kogut et a/.1996, Leitch et a/.1997, de 
Oliveira-Costa et aL1997, Jaffe et a/. 1998). Note, how- 
ever, that we assume no cross-correlation between free- 
free and dust, because any correlation makes the fore- 
ground separation easier. The same shape is also taken 
for synchrotron radiation. 

We choose amplitudes and frequency dependences for 
the galactic foregrounds consistent with the Kogut et 
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al. (1996) analysis of DMR, DIRBE and Haslam maps. 
We take the antenna temperatures of the free-free and 
synchrotron to vary with power-law indices -2.9 and - 
2.16, respectively. For the dust we assume a ly^ emissivity 
dependence and a single component with T = 18K. 

Draine and Lazarian (1997) have proposed an alter- 
native explanation to the observed correlation between 
dust and 30 GHz to 90 GHz radiation. They propose 
that the rotational emission from spinning dust grains, 
greatly increases the emission in the 10 GHz to 100 GHz 
range above what one expects from the vibrational emis- 
sion. We have not included this component of dust emis- 
sion in our model. Instead, we include something worse - 
a component with spectral shape similar to the "anoma- 
lous" emission, but which has no correlations with the 
dust. Again, this is more pessimistic than the strong 
correlation expected in a realistic model. 
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Fig. 1. — The frequency-dependent rms antenna tem- 
perature, giv\/l{l + l)Cii/(27r) evaluated at / = 500 (top 
panel) and I — 1500 (lower panel) for standard cdm CMB 
(black), dipole and thermal emission dust (both red), 
free- free (green), synchrotron (blue), SZ (cyan), and ra- 
dio and FIR point sources (both magenta). 

3.2. Extragalactic 

Extragalactic contributions to the microwave sky in- 
clude inverse Compton scattering of CMB photons by 
hot gas in clusters (the thermal Sunyaev-Zcldovich (SZ) 
effect), the Far Infrared Background (FIRB) and radio 
point sources. 

Following Tegmark and Efstathiou, we model the con- 
tribution from a diffuse background of unresolved radio 



point sources as having an intensity, cx with 

a white noise angular power spectrum [Ci independent 
of t) . Deviations from white noise due to clustering are 
observed to be neglible at 1.5GHz (TE96; TofoUati et 
al. 1998, herafter "To98"). Below 200 GHz, we take 
a = but above 200 GHz we introduce a break to a = 0.7 
as suggested by To98. We adopt this break in our spirit 
of pessimism because, despite decreasing the brightness 
of this contaminant, it actually makes determination of 
the CMB more difficult. This is due to the fact that with 
the break, the spectral shape more closely resembles that 
of the CMB. 

We are actually considering the power spectrum of the 
sources which remain after some cut is done of clearly 
contaminated pixels, e.g., those above a ba threshhold 
where ct^ is the variance in the map. Thus the amplitude 
depends both on the number-count distribution and on 
the level of flux cut that is used. Although this flux 
cut will vary for maps at different frequencies and from 
different instruments, we choose to fix it at 1 Jy. We view 
this as quite conservative since the typical level for all the 
Planck maps is about cr = 0.1 Jy. This is according to 
Tegmark & de Oliveira-Costa (1998) who used the point- 
source model of To98 and included the effect of reduction 
in cr that one can achieve by applying a Wiener filter. The 
values of a for the MAP maps should not differ by more 
than a factor of 2 from those for the LFI. 

For the amplitude of the FIRB we rely on the esti- 
mates of BG98 which are derived from the model of FIR 
point sources of Guiderdoni et al. (1998). This model has 
succesfully predicted source counts in a wide wavelength 
range, from 15 to 850 microns (see BG98 and references 
there in). The mean frequency dependence of the model 
is shown in Fig. |l|. Bouchet and Gispert (1998) have 
shown that this frequency dependence has only slight 
spatial variations, lending credence to our modeling of 
it as a frequency dependence times a fixed brightness 
spatial template. We assume clustering is unimportant 
and therefore the spatial power spectrum has the same 
shape as we have assumed for radio point sources: Ci is 
a constant. 

CMB photons moving through a hot gas of electrons 
have their frequency shifted as they Compton scatter, 
leading to the generation of anisotropy with a non- 
thermal spectrum. This Sunyaev-Zeldovich (SZ) effect 
can also be treated as an additive foreground, with the 
frequency-dependence of a Compton y distortion. Calcu- 
lations of the power spectrum of this foreground compo- 
nent, assuming a Press-Schechter distribution of clusters 
with masses greater than some cutoff have been done 
(Aghanim et al. 1996, herafter A96;Atrio-Barandella & 
Mucket 1998). We use the results of A96 for the 17 = 1 



4 



cosmology. Their power spectrum is well-fit in the range 
of interest by C; = a(l + Ic/l) where Ic = 1190 and a 
is such that l{l + 1)Ci/{2tt) = S.S^K^ at I = 1500 in 
the Rayleigh- Jeans (low frequency) limit (see Fig. |^). 
Modelling of this contribution will soon be improved 
by replacement of the use of Press-Schechter with N- 
body/hydro simulations. 

3.3. Spectral Shape Uncertainty 

Implicit in our formalism is that the frequency depen- 
dence of the foregrounds is known perfectly and has no 
spatial variation. However, we can allow for some de- 
gree of spatial dependence of the spectrum as follows. 
Consider foreground i with mean power-law frequency 
dependence, /?, and deviation SPim- Then, the signal 
contribution to the data, A^im from component i is 

(11) 

Thus we can treat radiation from a component with spa- 
tially varying spectral index as due to two components 
with amplitudes aum and aumSPirm which will, in gen- 
eral, be correlated. For simplicity we have modeled these 
additional components as uncorrelated with the parent 
component and taken {aumSpimau'm,' Pi'm') = Cii{5(3'^). 
We have assumed {5(3"^) = 0.25 for the rotating small dust 
grains, dust, and synchrotron with the same prior as used 
on other foregrounds. TE96 also considered using extra 
components to model spatial dependence of the spectral 
shape. For an alternative approach, see Tegmark (1997). 

3.4. Foreground Polarization 

Precious little is known about the polarization of fore- 
grounds. For a review, see Keating et al. (1998). Extrap- 
olation from observations at low-frequency (< 1 GHz) are 
complicated by Faraday rotation along the line-of-sight, 
which is negligible at higher frequencies. Measurements 
at higher frequencies are in the galactic plane in dense 
star-forming regions (Hildebrand & Dragovan 1995) and 
are not expected to be representative of the statistics 
at high latitude. We make the same assumptions about 
foreground polarization as BPS. They neglect polariza- 
tion in all foregrounds except for synchrotron and dust. 
For the synchrotron, they take Cf' = 0.2Cf and for the 
dust they take the model of Prunet et al. (1998, here- 
after "PSB") (see also Sethi et al. (1998)). It must be 
kept in mind that the PSB calculation relies on indirect 
arguments and is therefore quite uncertain, as is the syn- 
chrotron model, as the authors readily admit. 

4. Application to Planned Experiments 



4.1. Temperature 

In Fig. ^ one can see that MAP's ability to measure 
the power spectrum is not significantly affected by the 
foregrounds below t ~ 500. Going to smaller values of 
H. we have greater frequency coverage, and greater ratio 
of signal to instrument noise. The only thing that gets 
slightly worse as £ decreases is the relative amplitude of 
the galactic foreground power spectra, but this effect is 
overwhelmed by the others. Of course going to higher 
I we have less frequency coverage and a smaller ratio 
of signal to instrument noise. The galactic foregrounds 
still do not become a problem though since their relative 
power continues to decrease. 

What does become a concern at higher I are fore- 
grounds with rising angular power spectra: radio point 
sources and the thermal Sunyaev-Zeldovich effect from 
galaxy clusters. These alone are responsible for the devi- 
ation of AC; from the no foreground case, visible in Fig. 

I 

The impact of the Sunyaev-Zeldovich component is 
worth exploring more. It is quite possible that the ac- 
tual amplitude is ten times larger than in our baseline 
model. The A96 calculation ignores the contribution 
from filaments — which may actually dominate the con- 
tribution from the clusters, and it ignores the clustering 
of the clusters. If we increase the power by a factor of 10, 
and relax the prior on it to a = 0.1 from a — 1.1, AC; 
doubles in the range from I = 400 to / = 700. On the 
other hand, if we increase the power by a factor of 10, 
and do not relax the prior, AC; only increases by a few 
per cent. What we learn from this is that having some 
constraints on the power spectrum of the SZ component 
can be just as important as the actual amplitude. 

The usefulness of prior knowledge of the SZ C; is en- 
couraging. It suggests that the analysis of MAP data 
can profit significantly from accurate theoretical predic- 
tions of the statiscal properties of the SZ component. It 
also suggests that measurements of the SZ component 
in much smaller regions of the sky, which roughly con- 
strain the power spectrum, can be beneficial to the anal- 
ysis of the full-sky MAP data. Such analyses should 
be possible from combining MAP data with datasets 
from higher frequency instruments such as TopHat^ and 
BOOMERANG]], which by themselves will be extremely 
interesting CMB datasets. 

Planck's ability to measure the power spectrum is not 
significantly affected by the foregrounds below t ~ 1200. 



^ TopHat home page: ittp : //topweb.gsfc .nasa.gov 
\ ROOMER A NC home pagn: 



http : //astro . caltech . edu/~mc/boom/boom . html 
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At higher £, the frequency coverage reduces, the noise 
in each channel increases and the SZ, radio and FIRB 
components increase in amphtude. Unhke for MAP, SZ 
is not important because in the HFI frequency range, 
SZ is easily distinguished from CMB; there is even the 
null at 217 GHz. However, the radio point sources and 
FIRB are a concern. There is strong dependence on the 
prior. Even with moderate prior information (a = 1.1 
on these two components), ACq; is 3 times larger than 
the no foreground case. With an infinite prior this re- 
duces to a much less significant factor of about 1.2. The 
situation is greatly improved if the flux from the two 
backgrounds of unresolved sources is a factor of 4 less in 
amplitude (16 in Ci) than we have assumed. This is not 
unlikely since our assumed flux cut of 1 Jy is about 20 
times the level of confusing noise, calculated by Tegmark 
& de Oliveira-Costa (1998), in the (post- Wiener filter- 
ing) 143 GHz, 217 GHz and 353 GHz HFI maps, and is 
therefore an extremely conservative 20a cut. Thus, we 
also show the results with our input power spectrum for 
point sources, and the FIRB each reduced by a factor of 
16 as the dashed line in Fig. ^. 

We see that with only the use of a moderate amount 
of prior information, the errors on the Co;s here are 
not qualitatively different from the no-foreground results. 
The conclusions of those forecasting cosmological param- 
eter errors would not be qualitatively changed by includ- 
ing the effect of the foregrounds as modelled here. 

If galactic foregrounds are well-described by the model 
used here, then they will not have significant impact on 
the primary science goals of MAP and Planck. That is 
perhaps the most robust conclusion to draw from the 
above. This is not to say that these foregrounds do not 
have their impact on how well the CMB can be measured. 
The left side of Fig. shows how the foregrounds affect 
the uncertainties in a|^^. As long as 6afJf^/^/Ci < 1 
then sample- variance dominates the errors in C; . As can 
be seen in the figure, this inequality holds out to at least 
I = 500, except for MAP in the case of pixel-by-pixel 
subtraction (/? = 0, or no use of prior information). 

4.2. Polarization 

The CMB is expected to be polarized at a level of 
about 10% of the anisotropy. The polarization fore- 
grounds are no where near as well-understood and ex- 
plored as the temperature foregrounds. However, taking 
some initial guesses at the polarization foregrounds we 
find the outlook for CMB polarization measurement by 
MAP and Planck to be fairly bright. The reason being 
that, once again, there is a window in frequency space 
where the CMB is the dominant contributor to spatial 
variations in polarization. 
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Fig. 2. — MAP uncertainties for a = 0.1 on all fore- 
ground components (top curve), a — 0.1 on all fore- 
ground components except for radio point sources and 
SZ, for which a = 1.1 (second to top curve). The lowest 
uncertainty curve is identical for a = oo and no fore- 
grounds. Planck uncertainties for a = 0.1 on all fore- 
ground components except radio point sources and the 
FIRB for which a = 1.1 (highest curve); a = oo (middle 
solid curve) the no foreground case (bottom curve), and 
same as the top curve but with the FIRB and radio point 
source power spectra reduced by sixteen (dashed curve). 
With the FIRB and radio point source power spectra re- 
duced by a factor of 16, the a = oo case is identical with 
the no foreground case. 
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Fig. 3. — MAP (top) and Planck (bottom) component 
map uncertainties expressed as Saum/ VCa for the CMB 
(left) and the three galactic foregrounds (right). No other 
foregrounds were included in the calculation of the un- 
certainty. On the left, the three cases are pixel by pixel 
subtraction (/3 = 0, top curve), Wiener filtering (/3 = 1, 
middle curve) and no foregrounds (bottom curve). On 
the right the three cases are dust (red), free- free (green) 
and synchrotron (blue) with a. (3 — \ prior applied to all 
the components except for the CMB and the component 
in question. 



This window does not necessarily exist across the en- 
tire polarization power spectrum, and in particular may 
disappear at low /. This is unfortunate since the two 
most interesting features in the polarization power spec- 
tra are the bump at I ~ 2^Zri where the redshift of 
reionization and the B-mode power spectrum due to ten- 
sor and vector modes (Kamionkowski et a/.1997, Seljak 
& Zaldarriaga 1997) which also peaks at low I. 

Here we focus on the reionization bump. To study 
sensitivity to it we have not implemented Eq. ^. Instead 
we have ignored cross-correlations between temperature 
and polarization so that Eq. ^ is applicable with appro- 
priate substitutions (e.g., Cu — > C^f ). In general the 
cross-correlations improve the constraints on the polar- 
ization power spectrum (BPS) but that shouldn't be the 
case here since the reionization bump is a feature solely 
of the polarization maps and docs not show up in cross- 
correlation with temperature maps. 

For standard CDM with an optical depth to Thom- 
son scattering of t = 0.1, Planck measures the reion- 
ization feature with cosmic variance precision (although 
HEX alone does not and neither does MAP). At larger £, 
where the signal is large, the foregrounds, as modelled 
here, have no significant impact on the ability of either 
of the satellites to measure the CMB polarization power 
spectrum. Our infinite foreground prior (or Wiener filter) 
results are in agreement with the Wiener filter results of 
BPS. 




10 100 

1 



Fig. 4. — Cf vs. £ for standard cdm (thick black line) 
with optical depth to Thomson scattering of r — 0.1 and 
expected uncertainties under different assumptions. The 
assumptions are, from top to bottom in each panel, no 
foreground prior (green), infinite foreground prior (blue) 
and no foregrounds (red). 
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5. Discussion 

We have presented a method to calculate the sensi- 
tivity to the CMB and its power spectrum given multi- 
resolution, multi-wavelength observations of a sky that 
consists of multiple foreground contributions^. The ap- 
plications to MAP and Planck have allowed for much 
greater freedom in the behavior of the foregrounds than 
did previous analyses (TE96, BG98). Despite this extra 
freedom, the conclusions are similar — that foregrounds 
are not likely to qualitatively affect the uncertainties that 
can be achieved on cosmological parameters. Similar 
conclusions have been reached by de Oliveira-Costa et 
a/.(1998). 

Our approach has not fully taken into account the 
non-Gaussianity of the foregrounds, spatial dependence 
of the spectrum of each component, uncertainty in the 
spectral shapes, and unknown components (e.g., a popu- 
lation of points sources whose spectra peak at 90 GHz). 
For these reasons it is difficult to conclude with certainty 
that the foregrounds will not qualitatively affect the de- 
termination of cosmological parameters. However, a very 
important reason for our rosy conclusions is a very sim- 
ple fact: for most of the multipole moments measured 
by a given experiment, the quality of the CMB map can 
be highly degraded, without having any impact on the 
quality of the power spectrum. Thus, any effect we have 
not included here has to overcome this hurdle in order to 
be important. 

Non-gaussianity is both a friend and a foe. We have 
already exploited it here in assuming that the brightest 
points sources could be identified with threshhold cuts 
and removed. However, it can present a challenge to the 
above sample- variance argument if it resulted in the er- 
rors in each aim being strongly correlated with the errors 
in aim' , in such a way that they did not beat down with 
many averagings. One can think of this as an effective 
reduction in the number of independent modes in the 
foreground (Eisenstein 1998). However, we expect that 
small-scale behavior in patches of the sky sufficiently sep- 
arated to be decorrelated. Hence we do not expect the 
mode number reduction to be large, though further in- 
vestigation of effects of non-Gaussianity is clearly war- 
ranted. 

We have also negected things that will improve esti- 
mation of the CMB from MAP and Planck data, such as 
the use of maps at other frequencies, e.g., DIRBE, IRAS 
and FIRST (which will fly with Planck). Assumptions 
about the smoothness of foreground power spectra are 



IDL programs implementing this procedure are available from the 
author. 



also reasonable and could significantly reduce our error 
forecasts at high Z, by extending the information gained 
at lower I where there is greater frequency coverage. 

It is clear though, that even if foregrounds do not do 
anything more than double the errors on cosmologcial 
parameters, the determination of the exact size of the er- 
ror bars will probably be dominated by foreground con- 
siderations. Small patches of the sky will be analyzed 
separately, with those appearing the cleanest given more 
weight. Foreground model residuals will be agressively 
sought. Thus the study of foregrounds remains very im- 
portant. We close by listing the following improvements 
in our understanding of foregounds which could prove to 
be extremely beneficial: 

•More accurate theoretical calculation of the statistics of 
the SZ component. Our positive conclusions for MAP 
depend on the amplitude of the SZ power spectrum and 
on how well that power spectrum can be determined a 
priori. We have shown that having a prediction of Cf ^ 
good to about a factor of 2 (which would justify our use 
of a = 1.1) is enough to keep AC/ within about ten per 
cent of the no foreground case, even if Cf ^ is ten times 
larger than the A96 calculation. 

•Higher frequency complements to MAP, such as are 
coming from balloon flights (e.g., TopHat and 
BOOMERANG). Even coverage of just a few per cent of 
the sky, can be used to characterize the level of contam- 
ination in the rest of the sky. 

•A point source survey near 90 GHz (see Gawiser et 
al. 1998). 

•Further development of methods for removing non- 
Gaussian foregrounds and understanding of resulting 
CMB uncertainties. 

•A full-sky, high resolution Ha survey, since this is a 
tracer of free-free. Useful steps in this direction are al- 
ready being made with a survey in the North with one 
degree resolution (Reynolds et aL1998) and one in the 
South with 0.1 degree resolution (McCuUough et a/. 1998) 
nearing completion. 

•Measurements of high galactic latitude dust and syn- 
chrotron polarization. 

I thank K. Ganga, A. Jaffe and J. Ruhl for useful 
conversations, as well as the organizers and participants 
of the Sloan Summit on Microwave Foregrounds, who 
have informed the above discussion. I acknowledge the 
use of CMBFAST (Seljak & Zaldarriaga 1996). 
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